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ABSTRACT 

Context. This study is the third of a series that investigates the degeneracy and stochasticity problems present in the determination of 
physical parameters such as age, mass, extinction, and metallicity of partially resolved or unresolved star cluster populations situated 
in external galaxies when using broad-band photometry. 

Aims. This work tests the derivation of parameters of artificial star clusters using models with fixed and free metallicity for the 
WFC3-I-ACS photometric system. Then the method is applied to derive parameters of a sample of 203 star clusters in the Andromeda 
galaxy observed with the HST. 

Methods. Following Papers I & II, the star cluster parameters are derived using a large grid of stochastic models that are compared to 
the observed cluster broad-band integrated WFC3-I-ACS magnitudes. 

Results. We derive the age, mass, and extinction of the sample of M31 star clusters with one fixed metallicity in agreement with 
previous studies. Using artificial tests we demonstrate the ability of the WFC3-I-ACS photometric system to derive the metallicity of 
star clusters. We show that the metallicity derived using photometry of 36 massive M31 star clusters is in a good agreement with the 
metallicity previously derived using spectroscopy taken from literature. 

Key words, galaxies: individual: M31 - galaxies: star clusters: general 


1. Introduction 

Several recent studies (e.g., Fouesneau & Lan 9 on 2010[ Popescu 


& Hanson||2010[ [Anders et al.||2013[|de Meulenaer et al.||2013T 

hereafter Paper I) were aimed at the determination of partially 
resolved or unresolved star cluster parameters based on broad¬ 
band photometry and using more elaborated models than the 
traditional simple stellar population (SSP) models, strongly bi¬ 
ased by the presence of stochasticity problem in star clusters. 
The stochastic presence of massive stars in a star cluster causes 
its integrated colors to be dispersed, making SSP model predic¬ 
tions biased, especially in case of low-mass star clusters. Using a 
grid of star cluster models in which stellar masses are randomly 
sampled to take stochasticity into account, de Meulenaer et al. 
( |2014| hereafter Paper II) showed how it is possible to determine 
the parameters when the metallicity of clusters is unknown, and 
also quantified the accuracy of the metallicity derivation itself. 

In this study we apply these stochastic models to real clus¬ 
ters from the Andromeda galaxy. Using the Wide Field Camera 3 
(WFC3) and the Advanced Camera for Surveys (ACS) on board 
the Hubble Space Telescope (HST), the Panchromatic Hubble 
Andromeda Treasury (PHAT) team (see, e.g., jPalcanton et al.| 
2012[ Beerman et al.|2012 Weisz et al.|2013| l performed a pho 


tometric survey of 1/3 of the M31 galaxy, providing a large cat¬ 
alog of new clusters. [Fouesneau et al.j ( [2014| l derived the age, 
mass, and extinction of the PHAT “year 1” cluster catalog using 
the aperture photometry data provided by [Johnson et al.[ ( |2012| l. 

The objective followed here is to study a sample from the 
PHAT cluster catalog allowing variable metallicity in the model 
grid to explore the metallicity effects by comparing the results 


to the case with the fixed metallicity. We also aimed to derive 
the metallicity of objects with reasonably good photometric ac¬ 
curacy. 

The structure of the paper is the following: Sectionj^presents 
the method of star cluster parameter determination and tests its 
accuracy using the WFC3-I-ACS photometric system. Section 
[^presents the application of the method on the sample of the 
PHAT star clusters. Conclusions are presented in Section]^ 


2. The method of parameter derivation 

2.1. Presentation of the method 


As presented in Papers I & II, the determination of the physical 
parameters (age, mass, extinction, and metallicitjQ of a given 
observed star cluster is based on a comparison of its integrated 
broad-band photometry to a grid of star cluster models. 

The first task is thus to build a 4-dimensional grid of mod¬ 
els for every value of the four physical parameters. To take the 
stochasticity problem into account, we computed the grid of dis¬ 
crete cluster models by randomly sampling the stellar mass ac¬ 
cording to the initial mass function (IMF, Kroupa|200I i follow¬ 
ing the method described in Deveikis et al.| ( 2008[ ) (see also |San-| 
[tos & Frogel|1997|[Cervino et al.|2002| l. The stellar luminosities 
are extracted from isochrones of the selected age and me tallicity 
of the model. We used the PADOVA isochrone^ from Marigo 


' We refer to extinction and metallicity as E(B - V) and [M/H] here¬ 
after _ 

^ FAUUVA isochrones from “(JMU 2.6”: fittp://stev.oapd.inat.it/cmd” 
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Fig. 1. Age (panel a), mass (panel b), and extinction (panel c) derived for a sample of 10 000 artificial star clusters with true mass logjofM/Ma) = 4 
using WFC3+ACS photometry with cr = 0.05 mag for the Gaussian photometric errors for each filter. The true metallicity of the artificial clusters 
is specified on the right and the metallicity of the grid used to classify the clusters is specified on the top. 



not yet included. The grid was built according to the follow¬ 
ing nodes: from log[Q(f/yr) = 6.6 to 10.1 in steps of 0.05, from 
logjQ(M/Mo) = 2 to 7 in steps of 0.05, and for 13 metallicities: 
from [M/H] = -1-0.2 to -2.2 in steps of 0.2. This results in a 
grid of 71 values of age, 101 values of mass, with 1 000 models 
per node, hence ~1 x 10® models for each metallicity. To limit 
the number of models that need to be stored in computer mem¬ 
ory, the extinction was computed when the observed cluster was 
compared with the grid of models. It ranges from E(B - V) = 0 
to 1 in steps of 0.01, therefore 101 values for the extinction. We 


used the Milky Way standard extinction law from Cardelli et al. 
( [T989l l. 


In a similar way to Fouesneau & Langon| ( |201()) l, |Fouesneau 
|et al.| ( |2014[ ), and Paper II, we evaluated the likelihood of each 
node of the grid to represent the observed magnitudes of the clus¬ 
ter. Within each node, we first computed the likelihood of each 
of the 1 000 star cluster models by 


Tmodel — I I ,- 

V^O-y 


exp 


(magy.ob, - magy ^„dei)^ 


2(^2 


( 1 ) 


where / stands for one particular filter, magy for the observed 
and model magnitudes in that filter, and F for the total number 
of filters. For example, F - 6 for the WFC3-I-ACS photometric 
system we use in this study. Then the likelihood of the node of 
age f, mass M, extinction E{B - V), and metallicity [M/H] is the 
sum of the likelihoods of its models. 


2.2. 1-metallicity vs 1-metallicity: exploration of the 
metallicity effects 

We apply the method of star cluster parameter derivation 
using the WFC3H-ACS photometric system: UVIS1/F275W, 
UVIS1/F336W, ACS/F475W, ACS/F814W, IR/FllOW, and 
IR/F160W to characterize the ability of this photometric system 
to derive the physical parameters of clusters. 

We generated three samples of 10000 artificial star clusters 
with random age in the range log[Q(f/yr) = [6.6,10.2], mass 
fixed to log[g(M/Mo) = 4, with random reddening in the range 
E(B - y) = [0, 1]. Each of the three samples has a fixed metallic¬ 
ity, [M/H] = 0, -0.6, or -2. Gaussian photometric errors of cr = 
0.05 mag were added to the magnitudes of the artificial clusters 
for each filter. 

In this test, we derived the age, mass, and extinction of clus¬ 
ters using model grids of fixed metallicity ([M/H] = 0, -0.6, or 
-2), and the results are displayed in Fig.[T] In panel a), we see 
that using a metal poor model grid to analyze a set of metal rich 
clusters can result in overestimated derived ages by about a fac¬ 
tor 10. This is the well-known age-metallicity degeneracy (see, 
e.g., |Worthey|1994]|Bridzius et al.|2008[ ). In panel b) and c), the 
selection of grid model metallicity shows also significant effects 
on the derived mass and extinction. 

Comparing Fig. □ to the pure optical case (UBVRI) or opti¬ 
cal -I- NIR case {UBVRIJHK) shown in Paper II, the presence 
of the UV filter (here UVIS1/F275W) significantly helps narrow 
the scatter in derived age, mass, and extinction for main diagonal 
panels (where the metallicity of the used grid is the same as the 
true metallicity) as well as for out-of-diagonal panels (where the 
metallicities are different). In the latter case, the similar metal¬ 
licity effects are seen as in Paper II, producing offsets on the age 
and mass, and increasing the scatter in the extinction. 


N 

Lnode it, M, E{B - V), [M/H]) = 2 L™del,« , 

n=l 


( 2 ) 2.3. 1-metallicity vs 13-metallicities: taking metallicity effects 

into account 


where N - 1 000, the total number of models contained in 
the node. The procedure is repeated for each node of the 4— 
dimensional grid, and the observed star cluster is then classified 
with the parameters of the node, which maximizes the quantity 


In this test we derive age, mass, extinction, and metallicity of 
the same artificial clusters of fixed metallicity, but here we use 
a model grid containing 13 metallicities (from [M/H] = -1-0.2 to 
-2.0) to check the ability of the method to constrain the metal¬ 
licity. We studied the clusters allowing three different amounts 
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Fig. 2. Age, mass, extinction, and metallicity derived for a sample of 10000 artificial star clusters of true mass logjQfM/Mg) = 4 , with cr = 
0.03 mag of Gaussian photometric errors (block of panels a), cr = 0.05 mag (block b), and cr = 0.1 mag (block c) using WFC3+ACS broad-band 
photometry. The true metallicity of the artificial clusters is indicated on the top of each column of panels. In each panel, the numbers indicate the 
percentage of the 10000 clusters with derived parameter in the region between the dashed lines, which is centered on the true parameter value. 


of photometric errors to study their influence on the derivation 
of parameters; Gaussian photometric errors of cr = 0.03, 0.05, 
and 0.1 mag were added to each magnitude. 

Fig.j^shows the results; the panel b) block can be compared 
to the results shown in Fig. 6 of Paper 11, as the same amount of 
photometric error has been added to artificial cluster photometry. 
The presence of the ultraviolet F275W filter of the WFC3-I-ACS 
photometric system is important as it helps to significantly re¬ 
duce the biases in parameter derivation, including metallicity, 
compared to the UBVRI photometric system. 

Concerning photometric errors, in Fig.|^we see a much bet¬ 
ter agreement between derived and true parameters for low and 
medium added photometric errors (panel a block ; 0.03 mag of 
photometric etTors, panel b block : 0.05 mag) than in the case of 
higher photometric errors (panel c block : 0.1 mag). 

In Fig.|^ for each parameter, we define regions around the 
true value to quantify the accuracy of the method of parameter 
derivation for the different cases of photometric etTors tested. 
These regions, enclosed between dashed lines in the figure, are 
considered regions of correct parameter derivation for this quan¬ 
tification of method accuracy. For the age, the region is defined 
up to 0.5 dex around the true value (see the dashed lines), for the 
mass up to 0.2 dex, for the extinction up to 0.1 mag, and for the 
metallicity up to 0.5 dex. In each panel, we indicate the percent¬ 
age of the 10000 artificial clusters that have derived parameters 
in the correct regions. For example, in the case of photometric 
errors = 0.03 mag, 80% of the clusters with true solar metallicity 
have metallicity derived in the correct region; the percentages are 
60% for clusters with true [M/H] = -0.6 and 75% for clusters 
with true [M/H] = -2 (see bottom row of panel a block). When 
we increase the photometric errors to 0.05 mag (panel b block) or 


0.1 mag (panel c block), these likelihoods of correctly deriving 
parameters decrease strongly, especially for intermediate metal¬ 
licity. Indeed for the case of clusters with true [M/H] = -0.6, 
nearly half clusters have metallicities correctly derived when 
photometric errors = 0.05 mag, and only 30% when photometric 
errors = 0.1 mag. These numbers are only likelihoods marginal¬ 
ized for each parameter separately. 

We can also compute the joint age-metallicity likelihood, i.e., 
the likelihood of correctly deriving the age and metallicity pa¬ 
rameters simultaneously. For photometric errors = 0.03 mag, this 
likelihood is 80%, 60%, 75% for clusters with true [M/H] = 0, 
-0.6, -2, respectively. In this case of very high photometric ac¬ 
curacy, the age-metallicity likelihood is equal to the metallic¬ 
ity marginal likelihood (indicated in bottom of the panels block 
in Fig. 1^) because the age derivation is very accurate, and the 
limiting factor is thus metallicity. For photometric errors = 0.05 
mag, this likelihood is 75%, 45%, 65% for clusters with true 
[M/H] = 0, -0.6, -2, respectively. For photometric errors = 0.1 
mag, this likelihood falls to 65%, 25%, 50% for clusters with 
true [M/H] = 0, -0.6, -2, respectively. 

Hence, the better the photometric accuracy is, the higher the 
accuracy of the derived parameters can be achieved with the 
WFC3H-ACS system filters, especially for the metallicity, as will 
be emphasized in the following section with real data. 

3. Application to the MSI PHAT star clusters 

3.1. Cluster sample 

Using the WFC3-I-ACS photometric system on board HST, the 
Panchromatic Hubble Andromeda Treasury (PHAT) team (see, 
e.g., [Dalcanton et al.|[2012[ [Beerman et al.|p()12[ |Weisz et al.| 
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Fig. 3. Photometric accuracy of the sample of 402 clusters (all dots) selected with available photometry in all filters from the Johnson et al. j2012^ 
catalog. Group 1 clusters (65 objects, large filled dots) have upper limit of the accuracy indicated by long dashed lines in different filters. Group 
2 clusters (138 objects, large open dots) have accuracy <0.5 mag (short dashed lines) in all filters. Small gray dots designate clusters which are 
not members of group 1 or group 2 and not studied in the following because of too low photometric accuracy (they have accuracy > 0.5 mag in at 
least one filter). In each panel we indicate the number of clusters (from the total 203 clusters of group 1 and group 2) with photometric accuracy 
better than 0.1, 0.05 and 0.03 mag in the filter associated with the panel. 


|2013| l performed a survey of 1/3 of the M31 galaxy, providing 
a large catalog of 601 clusters pohnson et al. |2OT^. This cat- 
alog of star clusters was already analyzed by Fouesneau et al. 
( |2014| l, who derived their age, mass, and extinction. They used 
a constant solar metallicity through the whole M31 disk, argu¬ 
ing that the HII zone study (Zurita & Bresolin |2012[ ) does not 
show any significant metallicity gradient. They allowed four dif¬ 
ferent metallicities ([M/FI] = -0.7, -0.4, 0.0, and h- 0.4, i.e., 
from Small Magellanic Cloud to super-solar metallicities) for 30 
massive globular-like clusters with mass > 10 ^Mq. 


We first selected a sample of 402 clusters from the catalog 
of [Johnson et al. p012| l with available magnitudes in all pho¬ 
tometric filters. Then from this sample we selected two cluster 
groups, for which we display the photometric accuracy for each 
filter in Fig.[^ In cluster group 1 (65 objects, the large filled dots 
in Fig.[^, the photometric accuracy of objects is < 0.15 mag in 
F275W and F336W, < 0.1 mag in F475W and F814W, and < 0.2 
mag in FI lOW and F160W. In cluster group 2 (138 objects, the 
large open dots in Fig.[^, the photometric accuracy in each filter 
is <0.5 mag. In total, we analyze 203 clusters. In Fig. [^ we in¬ 
dicate for each filter the number of objects from the 203 clusters 
for which photometry is more accurate than 0.03, 0.05, 0.1 mag. 


We derived the parameters of both cluster groups 1 and 2 by 
firstly fixing the metallicity of all clusters to the solar value, and 
a second time by allowing a large range of metallicities in the 
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model grid with 13 values of [M/H] = -1-0.2 to -2.2, in steps of 
0.2. 

The color-color diagrams of cluster groups 1 and 2 are shown 
in Fig. 1^ in optical-only filters (panel a), and UV-optical-IR fil¬ 
ters (panel b), along with the SSP evolutionary tracks (tracing the 
age of clusters, valid for massive clusters only) of three metallici¬ 
ties as an illustration. In this figure, one can guess the advantage 
of using UV and IR photometry to derive the cluster metallic¬ 
ity. In panel a), with optical-only filters, we see that the SSPs are 
rather close, while this is not the case anymore in panel b), where 
UV and IR filters are shown. This indicates that, at least for mas¬ 
sive and old clusters (center to bottom of panel b), the derivation 
of metallicity is possible with the WFC3-I-ACS system, provided 
that the photometric accuracy is reasonable, such as for group 1 
clusters (for which the maximum photometric accuracy is indi¬ 
cated in the Fig.j^. 


3.2. Results with fixed solar metallicity 


The age, mass, and extinction derived with the fixed metallicity 
grid [M/H] = 0 are compared to the results of |Fouesneau et ah] 
( |2014| l in the first row of Fig.[^ A general satisfactory agree¬ 
ment is observed for the majority of clusters between the pa¬ 
rameters derived in this study and those supplied by [Fouesneau 


et al. (2014|l, except for a few objects. There are two reasons 


for the deviating clusters; (1) for 30 old and massive clusters, 
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Fig. 4. Color-color diagrams of the 65 group 1 clusters (filled dots) and 138 group 2 clusters (open dots) defined in Fig.|^ Panel a) shows the 
situation in optical colors only, while panel b) shows the situation when we make use of the ultraviolet and near-infrared filters. In both panels, the 
error bars indicate the maximum photometric error for the group 1 clusters. The three lines show SSP evolutionary tracks of metallicity [M/H] = 0 
(dotted line), -1 (dashed line), and -2 (solid line). The SSP ages extend from logjo(f/yr) = 6.6, in the upper part of each panel, to 10.1 in the 
lower part of each panel. The arrows indicate the direction of the extinction Ay = 1. 


[Fouesneau et al.| ( [2014| l derived parameters leaving the metallic¬ 
ity free to vary in the range [M/H]= [-0.7,-1-0.4], while we used 
fixed solar metallicity, and (2) the stellar models used to derive 
cluster parameters are very different. Indeed, PEGASE-basecj^ 
ste llar models ([ Eioc & Rocca-Volmerange|[T999[ l are based on 
the Bertelli et al.|(|1994|) stellar models w ith a simplified analytic 
description of |Groenewegen & de Jong| (jT99^ for the TP-AGB 
phase while PADOVA stellar models used here attempt to nu¬ 


merically reproduce the physics of that stellar phase (see Marigo 
|et al. poo's) [Girardi et al.poTOll. Eor the masses, [Fouesneau et al. 


( 2014| l provide the present-day mass, while the masses output 
by our method are initial masses, causing a natural slight shift 
between their mass predictions and our mass predictions. The 
four clusters of group 1 (Fig.[^; filled points) strongly below 
the ide ntity line and classified as massive by [Fouesneau et al. 


( 2014| l are in fact globular-like cluster^ wrongly classified by 
our method when the metallicity is fixed to the solar value in the 
model grid. We provide more details on them in the following. 


3.3. Results with free metallicity 

Parameters derived when metallicity is left free are shown in 
Fig.[^ As expected from artificial tests, the introduction of the 
free metallicity parameter introduces a new level of complexity. 
Here we compare the age, mass, and extinction derived when 
metallicity is left free versus the same parameters when metal¬ 
licity is fixed to the solar value. 

Many clusters seen as young or middle-aged in metal-fixed 
case are now seen as older. For example, the five points in the 
top left corner of the age panel (Fig.|^, four of which are the 
filled points below the identity line in Fig.[^) are seen as young 
in fixed solar metallicity case, but old in the free metallicity case. 
They are also more massive, located well above the identity line 
in the mass panel Fig.j^b). These five objects are classified as 

^ http://www2.iap.fr/pegase/ 

Verified inspecting the images provided by PHAT: 
http://archive.stsci.edu/pub/hlsp/phat/ 


low-metallicity when the metallicity is left free, which is also 
found likely when inspecting the images, as they are brighter 
in F275W and F336W filters than other globular-like clusters 
classified with higher metallicity. We used the individual clus¬ 
ter pictures in the six WFC3-I-ACS filters and also Sloan Digital 
Sky Survey (SDSS), Two Micron All Sky Survey (2MASS), and 
Galaxy Evolution Explorer (GALEX) images available through 
ALADIhQ Sky Altas to confirm that these objects are really 
globular-like clusters, and not young low-mass clusters. These 
five objects, with given PHAT ID 1439,428, 680, 683, and 1396, 
are also confirmed as globular clusters in the Revised Bologna 
Catalog (2012, version 5, see also|Galleti et al.|2004|l with desig¬ 
nation B064D-NB6, B229-G282, B165-G218, B167-G212, and 
NB21-AU5, indicated from 1 to 5 respectively in the Pigs. SHI 
and [8] 

As an illustration, 2-dimensional marginalized likelihood 
maps are shown for one of these objects, ID 428 (indicated as 
“2” in the Pigs. ss and[^, in Fig.|^ The likelihood maps are 
given for the cluster classification using all 13 metallicities of the 
model grid. The parameters derived taking the maximum of like¬ 
lihood Lnode (see Eq. in the 4-dimensional model grid are in¬ 
dicated with the white points in the 2-dimensional marginalized 
likelihood maps. Additionally, the black dots indicate the param¬ 
eters obtained when the metallicity is fixed to the solar value in 
the model grid (3-dimensional grid), resulting in a wrongly clas¬ 
sified younger, low-mass, and much more extincted solution. 

We note that in Pig. [^a) dozen of the objects with ages, de¬ 
rived assuming fixed solar metallicity, around logio(f/yr) = 8 
have overestimated ages when derived with free metallicities. 
This is very likely because these objects are relatively faint, with 
poor photometric accuracy, and are contaminated by bright red 
background stars. However, careful analysis of the WPC3H-ACS 
object images and their likelihood maps (similar to those shown 
in Pig.|7]l allows us to resolve degeneracies in most of the cases. 

Recently, Caldwell et al.| ( 201 1[ ) produced the spectroscopic 
study of old star clusters of M31 galaxy, using Lick indices to 


http://aladin.u-strasbg.fr/ 
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Fig. 5. Age, mass, and extinction derived with fixed solar metallicity grid vs the results of |Fouesneau et al^ ( |2014^ . Filled dots are the 65 clusters 
of group 1 and open dots are the 138 clusters of group 2, which are specified in Fig.[^ Clusters marked 1-5 in this figure and Figs.|^and[^are the 
same objects. 


Oi 







Fig. 6. Age, mass, and extinction using free metallicity grid vs the results derived using a model grid of fixed solar metallicity. Filled dots are the 
65 clusters of group 1 and open dots are the 138 clusters of group 2, which are specified in Fig.[^ Clusters marked 1-5 in this figure and Figs.|^ 
andj^are the same objects. 


derive their age, mass, extinction, and metallicity. To check the 
reliability of our derived metallicity, we compare those of the 36 
clusters common between the Caldwell et al. ( 201 l|l sample and 
the clusters analyzed in this study in Fig. ^ As Caldwell et al. 


( |201 l| l fixed the age of most clusters to 14 Gyr, here we only 
compare the mass and metallicity of the clusters. Again, an over¬ 
all agreement is found between the parameters. The accuracy of 
our photometric metallicity derivation is linked to the mass, and 
thus very likely to the signal-to-noise of available photometry 
for each object, as the scatter seen in Fig.[^c) is increasing with 
decreasing star cluster mass. The five clusters studied above are 
also indicated in Fig.[^ where one can see that the metallicity 
derived using our method coincides well with the spectroscopic 
method of |Caldwell et al.|P011| l. 


In Fig. 1^ we show the histograms of derived metallicity of 
group 1 (thick line) and group 2 (thin line) clusters for young 
age clusters (log[Q(f/yr) < 9, panel a) and old age clusters 
(logio(f/yr) > 9, panel b). Most of the young clusters are classi¬ 
fied as metal-rich, while the old cluster metallicities are more 
dispersed. Note, that even for low photometric accuracy. Fig. 
9(a) still tells us that the clusters are more likely metal-rich than 
metal-poor. Also, it is interesting to compare Fig. |^a) to the ar¬ 
tificial tests of Fig. in the worst photometric accuracy (panel 
c block). Indeed, we see a strong similarity between the case 
with true [M/H] = 0 and the left panel of Fig. Therefore, fix- 
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[M/H] [M/H] 


Fig. 9. The derived metallicity of clusters with ages lower than 1 Gyr 
(panel a), and of clusters with ages higher than 1 Gyr (panel b) for group 
1 (thick line) and group 2 (thin line) clusters. 


ing the metallicity to a value as high as the solar value assumed 
by [Fouesneau et al.| ( |2014| l is likely a good choice for most of 
the young M31 clusters. Note that Caldwell et al. ( 2009[ ) recom¬ 
mends supersolar metallicities for young M31 clusters, although 
they did not derive individual young cluster metallicities. 
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logio(Vyr) E{B-V) logio(il//Mg) 

Fig. 7. Two-dimensional marginalized likelih ood Lnode (see Eq.[^ parameter maps derived with free metallicity model grid for the cluster ID 428 
(B229-G282 in Revised Bologna Catalog V5, |Galleti et al.|2004[ indicated as “2” in the Figs.[^[^ and[^ photometry taken from [Johnson et af] 
\2Q\2) catalog. The white points indicate the maximum Lnode in the 4-dimensional parameter space, while the black points show the solution when 
the metallicity is fixed to solar value. The white contour lines contain 68% (solid line), 95% (long-dashed line), and 99% (short-dashed line) of the 
marginalized likelihood. 



l®Slo(-^-V-^-^©)caldwell P®/H] Caldwell l®Slo(-^-^/-^-^©)caldwell 


Fig. 8. Our results vs those of |Caldwell et al.| ( |2011^ for the 36 common massive clusters, shown for the mass (panel a), the metallicity (panel b). 
Panel c) shows the difference of metallicity between this work and |Caldwell et al.] l |201 1^ values vs the cluster mass. For the panel c), error bars of 
[Caldwell et al.H2011^ and derived here are summed quadratically. Filled dots are the group 1 clusters and open dots are the group 2 clusters, as 


defined in Fig.|^ Clukers marked 1-5 in this figure and Figs.[^and[^are the same objects. 


4. Conclusions 

We demonstrated, with the use of artificial star clusters, that the 
derivation of the age, mass, and extinction is possible when us¬ 
ing the WFC3-I-ACS photometric system, allowing us to avoid 
the biases introduced when the metallicity is fixed. In Paper II, 
we showed that optical (UBVRI) or opticald-NIR (UBVRIJHK) 


photometry cannot achieve the parameter derivation accuracy 
reached by using the WFC3-I-ACS photometric system; the im¬ 
provement in the latter system is because of the presence of the 
UV filter UVIS1/F275W of WFC3. In the case when the pho¬ 
tometric errors are reduced up to 0.05 mag for each filter, an 
accurate derivation of the metallicity parameter is also possible. 
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Using a sample from the M31 PHAT “year 1” star cluster 
catalog of Johnson et al. p012|l, we derived their parameters us¬ 


ing fixed solar metallicity consistent with the Fouesneau et al. 


( 2014| l study, and also with free metallicity in the model grid. 
For the clusters in common with the |Caldwell et al.| ( 201 l| l spec¬ 
troscopic study, we compare the metallicity with their values and 
found and overall agreement. This demonstrates, with the use of 
real star clusters, that the WFC3-I-ACS photometric system is fit 
to evaluate the star cluster parameters when the metallicity is un¬ 
known, and evaluate the metallicity when the signal-to-noise of 
photometry is high enough. 


Acknowledgements. We are grateful to the anonymous referee who helped im¬ 
prove the paper. This research was funded by a grant (No. MIP-074/2013) from 
the Research Council of Lithuania. 


References 

Anders, P., Kotulla, R., de Grijs, R., & Wicker, J. 2013, ApJ, 778, 138 
Beerman, L. C., Johnson, L. C., Fouesneau, M., et al. 2012, ApJ, 760, 104 
Bertelli, G., Bressan, A., Chiosi, C., Fagotto, F, & Nasi, E. 1994, A&AS, 106, 
275 

Bressan, A., Marigo, R, Girardi, L., et al. 2012, MNRAS, 427, 127 
Bridzius, A., Narbutis, D., Stonkute, R., Deveikis, V., & Vansevicius, V. 2008, 
Baltic Astronomy, 17, 337 

Caldwell, N., Harding, R, Morrison, H., et al. 2009, AJ, 137, 94 
Caldwell, N., Schiavon, R., Morrison, H., Rose, J. A., & Harding, R 2011, AJ, 
141,61 

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 
Cervino, M., Valls-Gabaud, D., Luridiana, V., & Mas-Hesse, J. M. 2002, A&A, 
381,51 

Chen, Y., Girardi, L., Bressan, A., et al. 2014, ArXiv e-prints 
Dalcanton, J. J., Williams, B. R, Lang, D., et al. 2012, ApJS, 200, 18 
de Meulenaer, R, Narbutis, D., Mineikis, T., & Vansevicius, V. 2013, A&A, 550, 
A20 

de Meulenaer, R, Narbutis, D., Mineikis, T., & Vansevicius, V. 2014, A&A, 569, 
A4 

Deveikis, V., Narbutis, D., Stonkute, R., Bridzius, A., & Vansevicius, V. 2008, 
Baltic Astronomy, 17, 351 

Fioc, M. & Rocca-Volmerange, B. 1999, ArXiv Astrophysics e-prints 
Fouesneau, M., Johnson, L. C., Weisz, D. R., et al. 2014, ApJ, 786, 117 
Fouesneau, M. & Lan 9 on, A. 2010, A&A, 521, A22+ 

Galleti, S., Bellazzini, M., & Ferraro, F. R. 2004, A&A, 423, 925 
Girardi, L., Williams, B. F, Gilbert, K. M., et al. 2010, ApJ, 724, 1030 
Groenewegen, M. A. T. & de Jong, T. 1993, A&A, 267, 410 
Johnson, L. C., Seth, A. C., Dalcanton, J. J., et al. 2012, ApJ, 752, 95 
Kroupa, R 2001, MNRAS, 322, 231 

Marigo, R, Girardi, L., Bressan, A., et al. 2008, A&A, 482, 883 

Popescu, B. & Hanson, M. M. 2010, ApJ, 724, 296 

Santos, Jr., J. F. C. & Frogel, J. A. 1997, ApJ, 479, 764 

Weisz, D. R., Fouesneau, M., Hogg, D. W., et al. 2013, ApJ, 762, 123 

Worthey, G. 1994, ApJS, 95, 107 

Zurita, A. & Bresolin, F. 2012, MNRAS, 427, 1463 


Article number, page 8 of[^ 













